Stiff-Glass Approximation of Mode-Coupling Theory 
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The stiff-glass approximation of tfie mode-coupling theory of the glass-transition — as introduced 
in a recent paper by Gotze and Mayr for a discussion of phenomena resembling the Boson Peak and 
High-Frequency Sound observations made in real glass-formers — is examined in detail. It amounts 
to a neglection of two-phonon processes and thus provides a clear physical picture for the MCT 
solutions deep in the glass. In addition, the effect of additional symplifying approximations and the 
combination of these approximations is studied. 
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I. INTRODUCTION 

In a recent paper, Gotze and Mayr demonstrated 
that the mode-coupling theory for the glass transition 
(MCT) appHed to a model of the hard-sphere system 
(HSS) at relatively high packing fractions yields spectra 
for the density correlation functions, that display a broad 
asymmetric peak as well as a sound-like mode at rather 
high wave vectors. Similar observations have been made 
in real glass formers and are known as the Boson Peak 
(BP) and high-frequency sound (HPS), respectively. To 
investigate these phenomena in detail, a large number 
of experiments has been performed, using various tech- 
niques such as Raman spectroscopy, x-ray spectroscopy, 
neutron scattering, or molecular-dynamics simulations. 
For an account of the literature, the reader is referred to 
the papers cited in This work also contains a sum- 
mary of the features signifying both the HPS and the 
BP. 

The mentioned aspects of the full numerical MCT so- 
lutions can be understood on the basis of two additional 
approximations; one being the generalized hydrodynamic 
(GHD) approximation, the other being the so-called stiff- 
glass approximation (SGA). Performing these two ap- 
proximations, one arrives at an equation that, while it 
still has to be solved numerically, provides physical in- 
sight on the origin of the BP- and HFS-like features of 
the HSS model. Especially the SGA can be understood 
on physical grounds as a description of phonon scatter- 
ing from the frozen-in amorphous structure, neglecting 
two-phonon processes. In a further simplification was 
carried out. It amounts to dropping the remaining wave- 
vector dependence of the microscopic frequencies in the 
GHD-SGA formula. As it leads to an equation that has 
been discussed earlier as a "schematic" toy model of the 
MCT, we shall call it the schematic-model description in 
the following. One obtains an analytically solvable ex- 
pression that, despite its simplicity, gives a good descrip- 
tion of the memory kernel responsible for the BP-like 
spectra. 

The aim of this paper is to drop this last simplification 
and to consider the SGA equations directly. This way, all 
of the additional simplifications applied on top of MCT 
can be investigated independently. We shall also compare 



with the analytic solution of ^ to understand its range 
of applicability. The discussion is meant to enlighten the 
physical mechanism behind the MCT boson peak, since 
the SGA is connected with a clear-cut physical picture. 

The paper is organized as follows: Section II will 
present a summary of the basic equations and approx- 
imations outlined above. In Sec. Ill, the solutions of the 



SGA equations with the full q dependence taken into ac- 
count will be presented. Section IV examines the effect 
of the GHD and schematic approximations, and Sec. ^ 
offers some conclusions. 



II. EQUATIONS OF MOTION 

A. Mode-Coupling Theory 

The mode-coupling theory of the glass transition 
(MCT) describes the dynamics of a glass-forming liq- 
uid in terms of the density correlators $q(t) = {Q*q{t)Qq)- 
Here, = exp(z(7r*j)/\/]V are the density fluctua- 
tions of wave vector g, for a system of N particles at 
density p. (...) denotes canonical averaging. For amor- 
phous systems, ^q{t) depends on the wave- vector only 
through its modulus g = At short times, one has 



^,{t) = Sq 



2 ^q^qi 



. , where Sq = {\q^\^) is the 



static structure factor specifying the equilibrium struc- 
ture of the system. $7^ = v^q^ / Sq are the phonon fre- 
quencies. A system of hard spheres shall be considered. 
The units of time and length are chosen in accordance 
with Q such that the thermal velocity v = 2.5 and the 
sphere diameter d — 1. 

The normalized density correlators (j)q{t), given by 
^q{t) = Sq4>q{t), obcy the equations of motion obtained 
via a Mori-Zwanzig projection operator formalism, 



ni 



mq{t-t')(j)q{t')dt' = 0. (la) 



Here mq(t) is the memory kernel, a correlation function 
of fluctuating forces. If one introduces Fourier-Laplace 
transformed quantities with the convention 4>q{z) = 



2 



; exp{izt)(j)q{t)dt, Eq. (la) is equivalent to 



(lb) 



where the fluctuation dissipation theorem has been used 
to connect (/)q(z) to the dynamical susceptibility Xq^z): 
Xq{z)l^q — zipq^z) + 1, with Xq denoting the isothermal 
compressibility. 

Within MCT, the memory kernel mq{t) is written as 
the sum of a regular contribution and a so-called mode- 
coupling contribution, the latter describing the slow dy- 
namics important close to the glass transition. For de- 
t_ails, the reader is referred to the original publications 
and standard literature One gets for mq{t) a 



quadratic fimctional in , m.q'^'^ [t) = Tq [ 



,MCT 



p=q-k (27I-)' 



:V{qkp)Mt)Mt), (2a) 



with coupling coefficients given by the equilibrium struc- 
ture of the system, 

Viqkp) = pSqSkSp[iqk)ck + {qp)Cpf/{2q^) . (2b) 

Here Cq denotes the Ornstein-Zernike direct correlation 
function, connected to the static structure factor by 
'^q — 1/(1 ~ f?Cg). If one introduces a wave- vector 
grid and approximates the integral as a Riemann sum, 
the mode-coupling functional can be written as J-q[f] = 
X^fcp Vqkpfkfp, where k and p run over the set of all grid 
points. The coefficients Vqkp are trivially related to the 
V{qkp). 

For the discussion of the MCT dynamics in the boson- 
peak region, Gotze and Mayr |^ chose a model of the 
hard-sphere system, given by evaluating Sq in the Percus- 
Yevick approximation, introducing a grid of M = 300 
wave vectors q with a cutoff of q* = 40, and by neglect- 
ing the regular part of the memory kernel. The depen- 
dence of the solutions on these approximations will have 
to be investigated seperately. For the structure factor 
used, a cutoff dependence of the solutions cannot be ne- 
glected, as was already mentioned in [|l[, and one is lead 
to the conclusion that the model described here is one 
of a soft-core fluid rather than true hard spheres. Even 
more so, the influence of the regular term is unclear at 
present. Recent calculations indicate that, at least in the 
white-noise approximation based on a generalized En- 
skog theory, the regular damping might be so large as 
to render the BP part of the spectrum nonexistent for 
a true hard-sphere system The question, however, 
remains unanswered still for regular potentials. Never- 
theless, since the concern of this paper is a comparison 
with the results presented in ||l[, we shall choose all nu- 
merical parameters the same. Equations (^ and ^ are 
closed then, and since Sq for the hard-sphere system does 
not depend on temperature, the only control parameter 
is the density p, which shall be written as the packing 
fraction ip = 7rp/6. 



B. Stiff-Glass Approximation 

The starting point for a theoretical understanding of 
the MCT high-frequency dynamics is a reformulation of 
the equations of motion, Eqs. (|l|), within the glass state. 
One introduces new correlators (t>q (t) that only deal with 
the decay relative to the frozen structure given by the 
nonergodicity parameter 4>q(t) = fq + {I - fq)4>q(t). 
Equations (|]) are covariant in the sense that the same 
equations hold for the new variables, 

Xqiz)/X°q = [z' - nlil - zmiz))] , (3) 

provided one replaces the static susceptibility by Xq = 
X° (1 - the frequencies fl^ by 

Cll - - /,) , (4) 

and the memory kernel by 

m,(0=.FW[<^]=.FW[<^]+.Ff [^], 
^'Hf]^J2^qkfk, 



^^Hf]=Y.^qkpfkfp. 

kp 

with new coupling coefficients given by 

Vqk=2{l-fq)J2Vqkpfp{l-fk), 
P 

Vqkp ^ (l - fq)Vqkp{l ^ .fk){l - fp) ■ 



(5a) 

(5b) 

(5c) 



(6a) 



(6b) 



The new variables (j)q (t) and rhq (t) have the properties 
that their long-time limits vanish, thus the low-frequency 
limits of their Laplace transforms are regular. 

To study BP phenomena, the limit of high packing 
fraction is of interest. There the original coupling co- 
efficients become large, formally written Vqkp — 0(1/77) 
with some small parameter 77. The virtue of the covariant 
transformation now is that it implies 1 — fq = 0{ri), and 
thus 



Vqk = 0(77) : 



Vqkp = 0(77') 



(7) 



i.e. the new coupling coefficients become small. This sug- 
gests to drop the second-order contribution to the mem- 
ory kernel rhq{t), such defining the stiff-glass approxima- 
tion (SGA), cf. Eq. (14) of |l|: 



.W(z) = -i/ [z-ny 



z + 



q,iSq (z) 



m^^\z) 



k 



Vqk4>^k\z) 



(8a) 
(8b) 



Equations similar to Eq. (^ were derived in a different 
context by other authors. Let us, for example, rewrite 
these equations in a notation closer to the conventions of 



3 



field theory. The dynamical susceptibility can be rewrit- 
ten for positive frequencies z = w^, a; > 0, into a prop- 
agator G{q,Lu) = ~x'q\y/z)/{xgCl'^), depending on the 

bare dispersion relation ^iq) — ^q{i + J2 k ^qk), and on 
the complex self energy w), to give 



G{q,Lu) 
T.{q,uj) 



1 



uj - e{q) - 
dkV,kG{k,iu), 



(9a) 
(9b) 



where Vqk = ^^Vgh^l are positive coupling constants. 
This is equivalent to the equations discussed by Martfn- 
Mayor et al. and Grigera et al. Q in connection with 
Euclidean random matrix theory, cf. Eqs. (13) and (17) 
of the latter paper, where, however, the vertices Vqk are 
different. Equations (||) can be viewed as a special case 
of the Dyson equation 

Equations (||) still have to be solved numerically. Thus, 
Gotze and Mayr ||l[ proposed two additional simplifica- 
tions, the first of which is the Generalized Hydrodynamic 
(GHD) description. One replaces the memory kernel 

■mq^\z) by its {q 0) limit, K{z) = m^q2Q{z). If one is 
only interested in a description of sound modes at small 
q, one could also replace the frequencies by their leading- 
order contributions, Ctq = Vocq + 0{q^), with the high- 
frequency sound velocity Wqo = v/ •\/(l — fo)So. In the 
following, we will only discuss the GHD-SGA obtained 
by keeping the full (/-dependence of Qq, 



X 



GHD 



ny z^-mi-zK{z)) 



(10) 



A further simplification was obtained in Q by re- 
placing VLq with some g- independent averaged value £7. 
This was motivated by the observation that in the range 
around the first diffraction peak in the structure factor, 
a q range known to constitute the dominant contribu- 
tion to the mode-coupling kernel [|| , the variation of the 
renormalized values Vtq with q is suppressed relative to 
that of Q,q. As a result, one arrives at 



4){z) = -l/ z-9?/ z + n^K{z) 
k(z) = wi0(z) , 



(11a) 
(lib) 



with wi ^ j dk Vok- The set of Eqs. ( 11a lib ) constitutes 
a well-known schematic model of MCT discussed earlier 
[Q, H, |lO|, the so-called Fi model. It can be solved triv- 
ially for either (j){z) or K{z). With the notation z = ui+iO 
we get 



/(2^), 
(12) 



with uj± = Izby^ilTj", and CJ^ = w^/fi^. The "semi-ellipse" 
spectrum if"(a;) thus obtained was demonstrated in 0| 



to have striking similarities to the spectra of the full MCT 
solutions in the BP regime at = 0.6. 

Eq. dH) implies a "gap" in the kernel spectrum for 
< w < UJ-. The full solutions show a pseudo-gap for 
frequencies below the BP frequency resembling this find- 
ing. The smearing out of the gap due to two-phonon pro- 
cesses can be understood on the basis of the schematic 
model by including these processes in a perturbative 
manner, which amounts to replacing K(lS) in Eq. (11a) 
with K'^'^\bj) — K{uj) + K'^'^^uj), where the latter term 
reads [01 



1 W2 



K{lu — uj')K" {ui') duj' . 



(13) 



with the integrated second-order coupling coefficient 
W2 = J dkVokk- K^'^\uj) is then stiU given by Eq. (p^), 
but one has to set (D^ = up' -f X^^^ (w). Analytic solv- 
ability can be regained by evaluating this in the (w — *■ 0) 
limit, ^ 0) = ir, 



W2 



w\ 37rrj 



(14) 



Since arriving at Eq. (|lj) involved three subsequent 
approximations which, in principle, could be applied in- 
dependently from one another, it seems worthwile to dis- 
cuss the effects of each approximation by itself. Espe- 
cially going over to the schematic model description of 
the memory kernel as a last step is, in contrast to the 
GHD-SGA, uncontrolled in the sense that its derivation 
is not based on a well-defined mathematical procedure. 

The GHD approximation was shown in [|l| to be of very 
good quality, provided on does apply it to the trans- 
formed memory kernel rhq^z). In other words, keeping 
the full (j-dependence of the nonergodicity parameters fq 
is of vital importance for a description of the BP and 
HPS spectra. It can thus be expected to be of similar 
quality when applied together with the SGA. 

While not leading to analytically solvable expressions, 
the SGA alone, Eqs. (^, has a clear-cut physical mean- 
ing: The first-order contribution to the memory kernel 
can be interpreted as an elastic scattering of phonons due 
to the frozen disorder in the glass, a mechanism thought 
to be dominant over two-phonon decay processes, which 
are represented by the second-order term and therefore 
ignored in the SGA. It should be noted, however, that 
this does not correspond to a harmonic approximation. 



III. SOLUTION OF THE SGA EQUATIONS 

To obtain numerical solutions of Eqs. (^), an algo- 
rithm adapted from the one used in ^ was employed. It 
solves the initial-value problem for the coupled integro- 
differential equations corresponding to Eqs. (|^) in the 
time domain. Regarding the stability of this method, 
proofs are available pTl], which apply to the SGA as a 
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FIG. 1: Transformed memory-kernel spectra rhq{Lo) of the 
fuU-MCT solutions for the HSS model, Eq. (||), for wave vec- 
tors q = 3.4, and q — 7.0, at packing fraction ip — 0.6 (solid 
lines). The dashed lines show the corresponding spectra ob- 
tained within the stiff-glass approximation (SGA), Eq. (H). 



special case. Here it merely is important that the mode- 
coupling vertices Vqk are positive. The solutions were 
Fourier transformed to give the spectra (t>q{uj) and m^'(a;). 

Let us start the discussion by a comparison of the mem- 
ory kernel spectra m'^{Lu), obtained from the solution of 
the fuU MCT equations, Eq. (|), and of the SGA equa- 
tions, Eq. (H). Figure Q shows such a comparison for 
packing fraction ip — 0.6. One notes a dominant broad 
and asymmetric peak responsible for the BP features of 
the correlator spectra The SGA solutions still re- 
cover the shape of that peak reasonably well; an obser- 
vation together with its form being reminiscent of the 
semi-ellipse known from the schematic model further jus- 
tifying the analogy put forward by the schematic model 
description. At frequencies u ^ 40, a pseudo-gap can be 
seen. As discussed in [Q, the SGA spectrum is not ex- 
actly zero within this gap; rather one obtains Rayleigh's 
law, TOq^-'"((jj) = RgLu^ + ©(w"'). The prefactor of the 
Rayleigh contribution is, however, so small that it is in- 
visible in the plot shown. Comparing with the full MCT 
solution, one clearly notes the filling of the gap due to 
the two-phonon modes. 

On the high-frequency wing, the spectra obtained from 
the full solutions show a nontrivial decay to zero, includ- 
ing, for larger q vectors, a second small peak at frequen- 
cies 200 ^ oj ^ 250. Interestingly enough, this peak is 
also reproduced by the SGA memory kernels. 

Since the SGA reproduces the memory kernel of the 
full MCT solutions resonably well in the frequency range 
of importance for the discussion of BP-like phenomena, it 
comes as no surprise that it also gives a basically correct 
description of the density correlator spectra ^'^{uj). Fig- 
ure ^ compares these spectra to those of the full MCT 
equations. The BP spectra are, with some quantiative 
differences, largely correct within the SGA, and the same 
holds for the HFS modes for those q vectors, where the 
sound damping is dominated by the BP background. At 
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FIG. 2: Dynamical structure factor S{q,u!) = Sq(f)'q{uj) for 
different wave vector q as labeled at packing fraction tp = 
0.6. The solid lines show the solutions obtained from the full 
equations, Eqs. (|^). The dashed lines are the corresponding 
SGA spectra, Eqs. (||). 

q ~ gD, where ~ 4 is the Debye frequency the SGA 
sound-mode damping clearly is too small. This failure is 
to be expected from the above discussion of the memory 
kernel, since in this case the sound resonance maximum 
occurs at frequencies above the BP spectrum, where the 
SGA memory kernel spectrum is too small. Similarly, 
the pseudo-gap at low frequencies and the rapid decay at 
large lo within the SGA are a result of the same features 
of the memory kernel spectra. 

The SGA is based on the two-mode contributions of 
order 0{rf') being neglegible in comparison to one- mode 
contributions of order 0{r]). The quality of this approxi- 
mation is expected to decrease with lowering the packing 
fraction, since rj increases with decreasing (p. Figures 
and ^ repeat the comparison made in Fig. |^, now for 
packing fractions ip — 0.5676 and ip — 0.54, respectively. 
One can identify the increasing intensity for low frequen- 
cies, eventually becoming so strong as to dominate over 
the BP background for Lp = 0.54. Nevertheless, the SGA 
gives the qualitatively correct picture even for (/? = 0.54 
at frequencies w 25, the number depending somewhat 
on the wave vector. One also notices the shift of both 
the sound resonance and the BP maximum to lower fre- 
quencies as the packing fraction decreases. This further 
demonstrates two features typical for the BP and the 
HFS 0. 
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FIG. 3: Comparison of full MCT and SGA for spectra ^g{Lo) 
as in Fig. ^, but for = 0.5676, and wave vectors q = 1.0, 
2.2, and 3.4. Other q were left out to avoid overcrowding the 
figure. 
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FIG. 4: Comparison of full MCT and SGA for spectra $'^{uj) 
as in Fig. ^, but for ip = 0.54, and wave vectors q = 1.0 and 
q = 3.4. 



One notes that for ip = 0.54 and (p = 0.5676 the true 
HFS maximum positions are higher than those read off 
the SGA solutions, especially for q vectors where the SGA 
significantly underestimates the damping. This effect is 
also known as level repulsion. For a discussion of the 
effects of the SGA on the HFS modes, Figs. and ^ 
show the maximum positions of the spectra (p'qluj) and 
the resonance widths versus wave vector q for different 
packing fractions. Two maximum positions are noted 
whenever the BP background and the HFS mode could be 
identified as distinct, separated by a spectral minimum. 
The full widths at half maximum (FWHM) were obtained 
by marking the lowest and the highest frequency where 
the correlator spectrum reaches half the value of its global 
maximum. 

For the HFS maxima as well as the BP maxima at 
ip = 0.6, Fig. 1^, good agreement between the full so- 
lutions, reproduced from and the SGA solutions is 
found. This even holds in the low-g regime, where only 
a damped sound mode can be seen in the spectra. Here 




FIG. 5: Maximum positions of the density-correlator spectra 
(j>'q{uj) obtained from the full MCT solutions at (p — 0.6 ver- 
sus wave vector q (open diamonds). Two symbols are drawn 
where two distinct maxima (separated by a minimum) could 
be identified. Filled circles indicate the maximum positions 
obtained within the SGA, Eq. (^), while crosses correspond 
to the results of the SGA-GHD, Eq. (^. The bars indicate 
the full widths at half maximum (FWHM) as described in 
the text. The smaller bars are FWHM of the SGA solutions, 
FWHM data for the SGA-GHD is left out to enhance read- 
ability. 



the shift in the maximum position due to damping is 
still small, such that an underestimation of the damping 
by the SGA does not necessarily lead to a wrong de- 
scription of the position. The agreement of the FWHM 
values also is good, with the reservation that the SGA 
values are too small in general. The lower frequency for 
which the spectrum reaches half the maximum intensity 
is described better than the higher one within the SGA, 
since here the pseudo-gap marks the dominant rise in the 
spectrum. 

Comparing the maximum positions and FWHM at </? = 
0.5676, one notes that for the BP maximum, the agree- 
ment between full MCT and SGA still is very good. Ex- 
pectedly, the description of the FWHM is worse, mostly 
due to the appearance of the pseudo-gap and due to the 
underestimation of the sound- mode damping around • 

The sound damping in the regime of the HFS is illu- 
minated in more detail in Fig. |7[ Here, full width at half 
maximum values F^ obtained from the full solutions of 
Eq. (H) and from within the SGA, Eq. (^), arc compared 
for the two packing fractions ip = 0.6 and ipi — 0.5676. 
For q < 0.6 (0.467) at ^ = 0.6 (0.5676), the SGA line 
widths could not be determined reliably. This reflects a 
limitation of the numerical procedure used, which is not 
adequate for very weakly damped oscillations, as occur- 
ing within the SGA. Nevertheless, the errors introduced 
by the numerical algorithm for low q do not significantly 
affect the accuracy of the solutions for larger q, since the 
small-(j contribution to the memory kernel is neglegible 

i- 

The dashed lines in Fig. u\ show the low-g asymp- 
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FIG. 6; Same as Fig. H, but for packing fraction ip = 0.5676. 
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FIG. 7: Double-logarithmic plot of linewidths (diamonds) 
and maximum positions $1™^'' (squares) of the density- 
correlator spectra versus wave vector q in the high-frequency- 
sound regime. The lower set of data corresponds to packing 
fraction — 0.6 (labels on left axis). The upper set shows 
the same data for (p = 0.5676 (labels on right axis). Open 
symbols are values determined from the solutions of the full 
MCT, filled symbols demonstrate the SGA solutions. The 
dashed lines correspond to the low-g asymptote, F^ = 5^7 
(see text). The dotted line corresponds to a power law, 
while the dot-dashed line demonstrates the Rayleigh contri- 
bution for ip — 0.6, multiplied by a factor of 10''. 



tote for the full MCT solutions, = 9^7, with 7 = 
''^oo'^5=o(^ ^ 0)- One notices an enhancement above 
the asymptotic law for q > 0.6. At around q — 1.4, 
the resonance width becomes of the same order of mag- 
nitude as its position, fl^'^^, and thus a lofFe-Regel limit 
is reached. Above this point, the line shapes differ signif- 
icantly from a Lorentzian and hybridization with the BP 
modes becomes important. For tp = 0.5676, the maxi- 
mum positions Q^^'^ are smaller than for ip = 0.6, while 
the damping changes only weakly in the regime of the 
enhancement over the asymptotic law, a behavior also 
found in experiment [ p^ p^ . Only for very small q vec- 
tors, q ^ 0.6, one notices in the full MCT solutions the 
expected behavior of Fg. 



In the q range from q « 0.6 to <? « 1.4, the data for 
tp = 0.6 can be approximately described by some power 
law, with its exponent being close to 4. To demonstrate 
this, the dotted line was included in Fig. 0, represent- 
ing a q"^ law. While in some investigations, no indica- 
tion of such a q* law was found [Q, it has been argued 
that this is behaviour is connected to Rayleigh's law . 
In the present case, however, the approximate q'^ power 
law is not an indication for Rayleigh scattering. As al- 
ready noted in ||l| , the Rayleigh contribution Rq is much 
lower. It is included in Fig. |^ as the dot-dashed line. Note 
that a magnification by a factor of 1000 was necessary to 
make this contribution visible in the plot. Thus one con- 
cludes, that the q'^ law for the high-frequency sound is a 
mere numerical coincidence. The possible appearance of 
a pseudo-power law for the HFS line width can be seen 
even more clearly for ip = 0.5676, where an exponent 
between 2 and 4 would be obtained. 

Fig. 1^ again demonstrates the validity of the SGA. At 
and below q ~ 0.6, deviations are to be expected, since 
the SGA line widths do not obey the asymptotic law. 
Rather one expects the SGA damping for 5 ^ to be 
given by Rayleigh's law noted above. Despite this fact, 
the overall description of F^ by the SGA is quite good. 



IV. EFFECT OF ADDITIONAL 
APPROXIMATIONS 

Having discussed the quality of the SGA in comparison 
to a full MCT solution, we shall now turn to a discus- 
sion of the additional simplifications outlined in the in- 
troduction. Since the q dependence of the memory kernel 
spectra is weak, cf. Fig. |^, one is lead to the combination 
of both the GHD approximation and the SGA. Figure ^ 
shows a comparison of thus obtained density-correlator 
spectra with selected full-MCT and SGA spectra repro- 
duced from Fig. |2[ Indeed, good agreement is found for 
the experimentally relevant q- vector range at q ^ ^d- 

The GHD-SGA solution fails to produce reasonable 
sound modes for wave vectors where the position of the 
sound maximum is at about the upper cutoff for the BP 
memory kernel. This is due to the sensitivity of the 
shape of the memory kernel in this frequency region. As 
noted in connection with Fig. |l|, the SGA still reproduces 
the nontrivial second peak in mq{uj) also seen in the full 
MCT solutions. Since this peak only evolves for higher 
q vectors, the GHD approximation obviously misses this 
feature, resulting in only weakly damped sound modes 
there. Still, the GHD-SGA combination is valuable for a 
discussion of BP-like spectra and the hybridization effects 
occuring between HFS modes and the BP spectra. Note 
that the maximum positions obtained from the GHD- 
SGA description show a qualitatively correct behavior. 
This can be seen from Figs. |] and |[ Within the hy- 
bridization regime, the FWHM values (not shown in the 
figures) also are almost quantitatively correct. 

A comparison between the full solutions and those ob- 
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FIG. 8: Spectra for <^ — 0.6 and q vectors as labeled ob- 
tained within the SGA (solid lines, reproduced from Fig. 0). 
Dashed lines show the corresponding spectra obtained within 
the GHD description in addition to the SGA. The dotted lines 
show examples for the full solutions (taken from Fig. 0). 



tained by applying the GHD only, still including two- 
phonon modes, was given in Fig. 6 of §. There it can 
be seen that the GHD approximation is, for q Qd, at 
least as good as the SGA, at the expense of still requir- 
ing the computation of a nonlinear memory kernel. Only 
the combination of both SGA and GHD seems to be too 
crude to deal with high-frequency sound modes for all 
wave vectors. For higher q vectors, the GHD-only ap- 
proximation of 0] and consequently also the GHD-SGA 
description somewhat miss the shape of the BP back- 
ground spectra, as is evident from the lower panel of 
Fig. ^. On the other hand, the SGA discussed here re- 
gains better quality for those ^'-vector regions where the 
HFS mode again shows hybridization with the BP-like 
spectrum. This shows that in this region, one-mode con- 
tributions from higher q vectors are important to fully 
explain the shape of the BP phenomena. 

As a last step, one can also compare the memory kernel 
spectra resulting from the SGA in the q ^ limit with 
the schematic model description proposed by Gotze and 
Mayr, Eqs. (0). Fig. 13 of [| presented a comparison 
at = 0.6 between the schematic model semi-ellipse and 
the fuU-MCT kernel spectrum at g = 0. Let us extend the 
discussion by adding the corresponding SGA spectrum. 
This is done in Fig. ||, while Figs. |l^ and ^ carry the 
same comparison to lower packing fractions, (p = 0.5676 
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FIG. 9: Memory kernel spectrum for (/? = 0.6 and q — 0. The 
spectrum m'^^o{u!) obtained from the full MCT equations, 
Eq. is shown as the thick solid line. The dashed line shows 
the SGA result. The schematic-model descriptions shown in 
Fig. 13 of |l| are also reproduced: the thin solid line is the 
semi-ellipse obtained from Eq. (^2|), the dotted line is the 
spectrum with r evaluated from Eq. (p^, see text for details. 
The dot-dashed line represents the result for the schematic 
model with a frequency-dependent second-order correction, 
Eq. (|l|), included. 
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FIG. 10: As Fig. |, but for ip = 0.5676. 

and (f> = 0.54, respectively. The averaged frequency was 
chosen to be 17 = 120 for (p = 0.6, in accordance with 
§, and = 85 (60) for (p = 0.5676 (0.54). One could 
improve on these values, but this would not give further 
physical insight, nor would it significantly alter the re- 
sults discussed below. 

The success of the semi-ellipse description at ip — 0.6 is 
based largely on the fact, that in the relevant q range, the 
variations of Clq are suppressed as rj becomes small. For 
lower packing fractions, the relative variations in fig are 
larger. Together with the increasing two-mode contribu- 
tions, this leads to stronger deviations of the fuU-MCT 
spectra from the schematic semi-ellipse form. 

As explained above, the two-mode contributions can 
within the schematic description be accounted for ap- 
proximately by including a frequency-independent damp- 
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FIG. 11: As Fig. a but for Lp = 0.54. 



ing term, Eq. (p^). While this procedure works well 
for ip = 0.6, it does not lead to an satisfactory result 
at the two lower packing fractions investigated here, cf. 
Figs. |l^ and |ll|. One needs to go beyond the white-noise 
approximation and include a frequency-dependent term, 
Eq. (|l^). Then the spectra are described fairly, but the 
advantage of analytic solvability is lost. 

However, the SGA solutions include no second-order 
contributions and thus should be compared to the pure 
semi-ellipse of the schematic description. It is remarkable 
that the SGA on the other hand catches the nontrivial 
shape of the BP-like contribution to the memory kernel 
spectrum qualitatively correct even at lower packing frac- 
tions, albeit with an overall too low intensity. The "tail" 
of the spectrum as well as the evolving central peak is 
missing in the SGA solution due to the well-understood 
occurence of the pseudo-gap. 



cussed and shown to provide a good description of the 
spectra in the boson-peak regime. It holds valid even 
for packing fractions lower than — 0.6, where the an- 
alytical formula of Gotze and Mayr Eq. (12), can- 
not describe the memory kernel any longer. This fail- 
ure of Eq. (12) is due to increasing variations with q 
in the phonon dispersion fig, which in the schematic 
model is neglected. The shape of the SGA memory ker- 
nel shows more structure than a simple semi-ellipse form. 
This leads to a substantially correct description of high- 
frequency sound and sound damping in this approxima- 
tion. 

Even the combination of the SGA with the general- 
ized hydrodynamic description (GHD) gives reasonable 
results for the memory kernels as well as for the density 
correlators, for q vectors where the damping of the sound 
mode is dominated by contributions from the BP-like 
memory kernel. Outside this window, one gets weakly 
damped sound waves, an artifact introduced in the in- 
termediate q range mainly by the GHD, not the SGA. 
Nevertheless, the GHD-SGA gives a basically correct de- 
scription of the BP spectra. This finding further corrob- 
orates the explanation of the HPS as a result of the BP 
spectra. 

As a side effect, it was demonstrated, that equations 
of the form of Eqs. (^) can be solved relatively easy in 
the time domain. At first glance, this seems paradox, 
since in Eq. (^), both q and frequency z only appear as 
parameters, while in the time domain one has to solve a 
convolution integral. 
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